Fractional quantum Hall states in lattices: Local models and physical implementation 
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The fractional quantum Hall (FQH) effect is one of the most striking phenomena in condensed 
matter physics. It is described by a simple Laughlin wavefunction and has been thoroughly studied 
both theoretically and experimentally. In lattice systems, however, such an effect has not been 
observed, there are few simple models displaying it, and only few mechanisms leading to it are 
known. Here we propose a new way of constructing lattice Hamiltonians with local interactions and 
FQH like ground states. In particular, we obtain a spin 1/2 model with a bosonic Laughlin like 
ground state, displaying a variety of topological features. We also demonstrate how such a model 
naturally emerges out of a Fermi-Hubbard like model at half filling, in which the kinetic energy part 
possesses bands with nonzero Chern number, and we show how this model can be implemented in 
an optical lattice setup with present technology. 



The FQH effect [1 is one of the most fundamental 
phenomena in strongly correlated electronic systems and 
a paradigm of topological behaviour. It has been thor- 
oughly studied and characterized in situations, where the 
solid surrounding the electrons has a modest effect on 
their properties [2H1]- This was possible, in part, thanks 
to the successful description in terms of the celebrated 
Laughlin wavefunctions [5] . The appearance of FQH be- 
haviour when the lattice structure created by the solid 
becomes important is, however, not so deeply studied nor 
understood. Nevertheless, the interest in lattice systems 
[51115] has been sparked by the possibility of observing 
such behaviour with ultracold atoms [MHTT] . obtaining 
high-temperature quantum Hall states [TS1424J . or even 
its applications in quantum information processing |25| . 

The standard route to the FQH effect in solids involves 
fractionally filled Landau levels (created by a strong, ex- 
ternal magnetic field) and electron-electron interactions. 
This mechanism can be adapted to lattice systems by 
replacing the fractionally filled Landau level by a suffi- 
ciently flat and fractionally filled Chern band and adding 
near- neighbour interactions [201 l^2H^ [251157] . A Chern 
band can be obtained by a proper choice of complex hop- 
ping amplitudes and can be made flat with longer range 
hoppings [IHUn]- 

Here, we propose and investigate a different mechanism 
to obtain FQH states in lattice systems. We consider 
spin 1/2 particles on a square lattice with some specific 
short-range interactions. As we will show, the ground 
state of the system is extremely well described by the 
Kalmeyer-Laughlin (KL) wavefunction [251 US] (the spin 
version of the bosonic Laughlin wavefunction with Lan- 
dau level filling factor i/ ~ 1/2) and exhibits the topologi- 
cal behaviour expected for FQH states. In particular, we 
consider periodic boundary conditions and analyze the 
ground state (quasi) degeneracy, its response to twisting 
the boundary conditions, the local indistinguishability of 
the ground states, and the topological entanglement en- 
tropy. We also demonstrate how our model can be im- 
plemented with ultracold atoms in optical lattices. The 
temperatures, tunnelling amplitudes, and interactions re- 
quired to observe the exotic topological behaviour are the 



same as the ones required to observe the Neel antiferro- 
magnetic order in the standard Fermi-Hubbard model, 
so that our predictions can be tested with present or 
planned technology [5DH55] . 

The mechanism leading to FQH states can be viewed 
from two perspectives. The first takes a spin wavefunc- 
tion with appropriate topological properties (like the KL 
state) that can be built out of correlators from a confor- 
mal field theory (CFT) (in our case the SU{2)i Wess- 
Zumino-Witten) and for which one can derive a (nonlo- 
cal) parent Hamiltonian. This Hamiltonian is then de- 
formed into a local one without crossing a phase transi- 
tion. The second perspective considers spin 1/2 fermions 
moving on a lattice with (not necessarily flat) Chern 
bands and strong local interactions. When the lattice 
filling is 1/2 for both spin up and down, a Mott state is 
formed, and the spin state inherits the topological char- 
acter of the Chern bands. The first view may be extended 
[T21[33] to other models to obtain, e.g., Moore-Read like 
states [3S]. The second lends itself for physical imple- 
mentations. 



I. MODEL 

We start out with a spin wave function for an even 
number N of spin 1/2 sitting at fixed positions in the 
two-dimensional plane 



N 



(1) 



n=l 



Here, z„ = a;„+zy„ is the position (a;„, y„) of the nth spin 
written as a complex number, s„ — ±1 labels the possible 
states ('up' and 'down') of the nth spin, and (5s = 1 for 
'J2n ■Sn = and 6s = otherwise. (IT]) is closely related 
to the KL state and reduces exactly to it in the case of 
a square lattice of infinite extent JJ^ . We shall refer to 
^ as the 'CFT state' because it can be expressed as a 
CFT correlator (see the supplementary information). 
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FIG. 1. Model and implementation, a, The considered spin lattice Hamiltonian is a sum of local two- and three-body inter- 
actions, b, For suitable parameters, the spin lattice Hamiltonian is effectively equivalent to the Fermi-Hubbard like Hamiltonian 
in equation S, of which we here show the kinetic energy parts Hkin.o-, o" G {tii}- Specifically, each arrow/line/wiggle from 
position n to TTi on the lattice represents the contribution t{e^'^ a]^„an<j -t-e^^^ajjo-amo-) to iifkin.o- with i and 4> given in the figure. 

c, A'^ fermions trapped in the optical lattice potential we propose to use for implementing the Fermi-Hubbard like Hamiltonian. 

d, We encode the spin up and down states in four internal hyperfine levels of the fermions. The blue/red states feel the blue/red 
potential in c and are hence trapped at the blue/red lattice sites. In this setting, we can implement the nearest-neighbour 
hopping terms through Raman transitions as indicated with the dashed blue and red lines in d. For this we need the three 
standing wave laser fields shown in b with electric fields: Er\ — Ke{—2Eez sin(fca;)e 
Er2 = Ke{—2iEez sm(ky)e~^"^'^) (the dashed red line along the y-axis), and Ebs ■ 

blue round shadow). Here, E = E* and E± are amplitudes, e^ = (0, 0, 1), £± = {^1, —i, 0)/\/2, k = uJa/c = n/a, where c is the 
speed of light and a is the lattice constant, fcz = uJt/c, ujr,b are frequencies, and T is time. The next-nearest neighbour hopping 
terms are implemented as a combination of hops induced by the fields listed above and tunnelling between nearest neighbour 
sites in the blue/red lattice. 



^'""^) (the dashed red line along the a;-axis), 
Re{2{E+e+ + E^e-)cos{kzz)e-^'^''^) (the 



An exact parent Hamiltonian of ippQ^^ has been de- 
rived in [12] by use of a quite general technique relying 
on CFT. The Hamiltonian consists of interactions be- 
tween all pairs and triples of spins in the system, which 
is challenging to achieve experimentally. One may argue, 
however, that deforming the Hamiltonian will not dra- 
matically change the physical properties of the ground 
state as long as no phase transition is crossed. With this 
in mind and specializing to an L^ x Ly square lattice 
with Lx even, we investigate the local Hamiltonian that 
remains after removing all long-range interactions and 
making all coupling strengths position independent: 



H — J2 2_^ 2S'„ • Sr, 
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Here, Sn = {3^,3^,8^) is the spin operator of the nth 
spin. As indicated in figure fli, the first (second) two- 
body term is summed over all pairs of nearest (next- 
nearest) neighbour spins, and the three-body term, which 
breaks time reversal symmetry, is summed over all trian- 
gles of neighbouring spins (for each triangle only one term 
is included and n, m, p label the vertices of the triangle in 
the counter clockwise direction as indicated with the ar- 
row). In this and the following section, we choose J2 = 1, 
J2 = 0, and J3 = 1/2 to make H as local as possible. The 
Hamiltonian conserves the total spin 3tot = X^n ^n and 



N La: X Ly 

12 4x3 

16 4x4 

20 4x5 
30 6x5 



Plane 

0.9860 
0.9812 
0.9728 



Torus 

K^Toi^?r)i KV'Tiii/'?r) 



0.9818 
0.9747 
0.9655 
0.9258 



0.9533 
0.9572 
0.9200 
0.9361 



TABLE I. Overlap between the ground state(s) of H 
and the CFT state(s). Columns 1-2: We consider A^ spins 
on a square lattice of size Lx x Ly. Column 3: The over- 
lap between the ground state -0po of the local Hamiltonian 
in equation (pi with open boundary conditions and the CFT 
state in equation (IIJ. Columns 4-5: On the torus, i.e. for 
periodic boundary conditions, there are two CFT states (see 
supplementary equations (S9-S10)), and the columns show 
the overlap of these states with the two lowest energy eigen- 
states ipTO and ipTi of S. The overlaps are remarkably high, 
in particular when taking the large dimension of the involved 
Hilbert spaces into account (e.g., 1.1 x lO" for a 6 x 5 lattice). 
We include fewer results for the plane than for the torus be- 
cause the lack of translational invariance in the plane makes 
it harder to diagonalize the Hamiltonian in this geometry. 



is hence 3U{2) invariant. 

To test how much the deformation of the Hamiltonian 
affects the ground state, we compute the overlap between 
the ground states before and after the deformation in ta- 
ble IT] column 3. For many-body systems, the dimension 
of the Hilbert space grows exponentially with the number 



of spins in the system, and the hugeness of the Hilbert 
space may easily cause states in the same phase to have 
poor overlaps. The fact that we obtain high overlaps for 
states with a few tenth of spins, on the other hand, is a 
strong indication that the deformation of the Hamilto- 
nian docs not bring the system to a different phase. In 
fact, we may regard (fTl) as an analytical approximation 
to the ground state, which is extremely helpful to see the 
underlying physics and to do computations for systems 
that are too large for exact diagonalization. We have 
chosen the parameters of the Hamiltonian to get large 
overlaps, but we note that high overlaps are obtained in 
a broad parameter region as we show below. 

The model we have looked at so far is defined in the 
plane and hence has open boundary conditions. To inves- 
tigate the topological properties of H and the CFT state, 
we shall also need to consider periodic boundary condi- 
tions, in which case H is translationally invariant. On 
the torus, i.e., for periodic boundary conditions in both 
directions, there are two CFT states, "0X0 """ ^"^"i V'Tf'^j 
for which we derive explicit analytical expressions in the 
supplementary information. Rewriting these expressions 
for the case of a square lattice, we find that the CFT 
states are exactly the KL states on the torus for all L^ 
and Ly. Table IT] shows that the large overlap between the 
ground states of H and the CFT states also hold on the 
torus. 



II. TOPOLOGICAL PROPERTIES 

FQH states are examples of topological states, and in 
the following we demonstrate that H and the CFT states 
exhibit topological properties that fit the properties of 
the V = 1/2 Laughlin state in the continuum. Charac- 
teristic features of topological states are a ground state 
degeneracy that depends on the topology of the surface 
on which the states are defined, the lack of ability to dis- 
tinguish the topologically degenerate states through local 
measurements, and a nonzero topological entanglement 
entropy. For FQH states in solids another important 
quantity is the Hall conductivity. The Hall conductiv- 
ity is closely related to the Chern number [55H55] , which 
can also be computed for spin systems [32] ■ For Laugh- 
lin states the expected behaviour is a many-body Chern 
number of the ground state manifold on the torus that 
is unity and spectral flow of the ground states into each 
other for certain choices of the lattice size under contin- 
uous twisting of the boundary conditions [231 HD] • 

We compute the spectrum on the torus through exact 
diagonalization [JT] of H (figure [2k) . The exponential 
increase of the dimension of the Hilbert space with the 
number of spins limits the system sizes we can consider. 
It is seen that the two lowest energy states approach each 
other as the system size is increased, and extrapolating 
the energy difference between the third lowest and the 
lowest energy state to the limit TV — )■ oo points to the 
existence of a gap in the thermodynamic limit. 
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FIG. 2. Topological properties of the Hamiltonian. a. 

Energy spectrum in the S^t = subspace on the torus for 
different lattice sizes (in green, A'^ = L^Ly) obtained by ex- 
act diagonalization. The number to the right of each energy 
level is the degeneracy and only the energies of the five lowest 
states {n — 0, 1, 2, 3, 4) are displayed. Note that the spectrum 
in the complete Hilbert space is the same except that each 
state is replaced by 25 -I- 1 degenerate states, where S is the 
total spin quantum number of the state. As the two lowest 
states have 5* = 0, the results suggest that there are two de- 
generate ground states and a gap to the first excited state in 
the thermodynamic limit like for the v = 1/2 Laughlin state 
in the continuum [47]. b. Energies of the five lowest states in 
the Slot ~ subspace on the torus for twisted boundary con- 
ditions in the a::-direction {Ox is the twist angle) and a lattice 
size of 6 X 5. Twisting the boundary conditions corresponds 
to gradually inserting a flux line through the hole of the torus 
(inset), and we observe that the two ground states flow into 
each other under this operation. 



For finite size systems, the local indistinguishability 
of topologically degenerate states is not perfect, but the 
ability to distinguish the states decreases exponentially 
with the size of the system. In figure |3j we plot the 
deviation d between the correlators of local operators 
computed for the two states on the torus for both the 
exact ground states and the CFT states. In both cases, 
d decreases with increasing system size, and for the CFT 
states the decay is clearly exponential. We also observe 
that, as far as the operators appearing in the Hamilto- 
nian are concerned, d is not small unless the lattice size 
is at least 6x5. This explains why the ground state de- 




FIG. 3. Indistinguishability of local observables. a, 

To demonstrate local indistinguishability of the states on the 
torus, we consider the set of all spin operators that act on a 
plaquette of four spins. Note that all plaquettes are equiv- 
alent due to the translational invariance. Using symmetries 
and the properties of spin operators, the correlators of all 
such local operators can be expressed in terms of the eight 
correlators depicted (the uppermost drawing, e.g., represents 
the correlator c{tp) = {ip\S^-^S^2\''P) / W''^) ^ where the spins 
on the plaquette are labelled ni, n2, n^, 714 in the counter 
clockwise direction starting from the upper right corner). For 
the special case Lx ~ Ly, the correlators displayed in black 
are not needed, b and c, Dependence of the relative dif- 
ference d — |2(co — ci)/(co + ci)| between the correlators 
Co = c{ipTo'^) and Ci = c{ip^f'^) on the size of the lattice 
for even-by-even and even-by-odd lattices (we use the mark- 
ers indicated in panel a). The extra set of smaller fainter 
symbols for 4 x 4 in b and 4 x 3, 4 x 5, and 6 x 5 in c (see 
the upper axes) show the same for the exact ground states of 
the Hamiltonian. The results are obtained by exact compu- 
tations for LxLy < 30 and from Monte Carlo simulations for 
LxLy > 36. The error bars of d are given as d ± Sd, where 

Sd = ^(||)2(&o)2 + (||)2(5ci)2, Sc, = Vvar(cO/AA, and 

the variance is taken over the outcome of A/" independent 
Monte Carlo trajectories. 
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FIG. 4. Topological entanglement entropy. To compute 
the topological entanglement entropy, we map the state in 
equation (II]) to a cylinder (inset) with closed ends at ±00 
(see supplementary information for details). We then cut the 

cylinder in two halves (part A and B) as indicated with the 

(2) 
blue plane, and compute the entanglement entropy S)^ = 

— ln(Tr(p^)) as a function of Ly for fixed Lx, where pA is the 
reduced density operator of part A. We choose here to use 
the Renyi entropy with index 2 because it is less demanding 

to compute numerically than the von Neumann entropy [48j . 

(2) 
For large Lx and Ly, Sj^ is independent of Lx and grows 

linearly with Ly, and the intersection with the vertical axis is 

the topological entanglement entropy (—7) |42l [^-52 . The 

fact that the results for Lx = 10 and Lx ~ 16 practically 

coincide shows that we are already in the Lx ^ 00 limit. The 

black dotted line is a linear fit to the results for Lx ~ 10 

and Ly > 10 and gives 7 = 0.374. The solid cyan line is the 

same fit except that 7 is fixed to the value 7 = ln(2)/2 « 

0.347 of the v = 1/2 Laughlin state, and it is seen that both 

fits fit the data well. The results for L^, = 4 are computed 

using the lowest energy state of the Hamiltonian in equation 

S with periodic boundary conditions in the y-direction, and 

they approximately follow the results for the CFT state. The 

results for Lx ~ 10 and Lx ~ 16 are computed from Monte 

Carlo simulations of r = exp{—Sj^ ). The error bars are given 



as S 



(2) 



±<\ 



where SSr 



j^^l Vvar(r)/AA and the 
variance is taken over the outcome of TV independent Monte 
Carlo trajectories. 



agreement with this value. For the exact ground state of 
H, we can again only consider rather small systems, but 
we observe that the results approximately follow those of 
the CFT state and also point to a nonzero value of 7. 



generacy first appears approximately for the 6x5 lattice 
in figure [2^, and the general decrease of d in figure |3] pro- 
vides further evidence that the two lowest energies will 
approach each other further for larger lattices. 

In figure |4] we compute the topological entanglement 
entropy —7 using the approach proposed in |42]. For the 
v = 1/2 Laughlin state in the continuum, 7 = ln(2)/2 
[33] , and the result we get for the CFT state is in perfect 



We have computed the many-body Chern number for 
both the exact ground states and the CFT states on the 
torus for a 4 x 5 lattice using the method described in (TS] 
155] . and in both cases we get 1. Figure [2t) demonstrates 
the flow of the ground states of H into each other when 
the boundary conditions are twisted, and for the CFT 
states this can be shown analytically. Again, our findings 
are thus in agreement with the properties of the t/ = 1/2 
Laughlin state. 



III. CONNECTION TO A FERMI-HUBBARD 
LIKE MODEL 

As mentioned in the introduction, each spin in our 
model may represent the spin of a fermion sitting on a site 
in an optical lattice. Let us consider the Fermi-Hubbard 
like Hamiltonian 



N 



HfH — 2_^ ^kin,ff + U 2_^ ^ii-f'^ 



"t'^H'^'4 



(3) 



on a square lattice, where ana annihilates a fermion with 
spin (J on lattice site n, U is a positive constant, -ffkin.o- 
is defined in figure [T|d, and we choose the total number 
of fermions to equal the number of lattice sites, ifkin.o- 
represents spin preserving hopping of fermions between 
nearest and next-nearest neighbouring sites with complex 
hopping amplitudes, and the second term in 7?fh repre- 
sents interactions between fermions sitting on the same 
site. 

When U is much larger than the hopping strengths 
t,i', we are in the Mott insulating regime, where each 
site is occupied by a single fermion. In this limit, we 
can use the Schrieffer- Wolff transformation [4¥ to derive 
an effective Hamiltonian Heg acting in the space of states 
with only single occupancy on all sites in the same way as 
the Heisenberg model is derived from the standard Fermi- 
Hubbard model. Reexpressing H^b in terms of spins, we 
get Hcs ~ H + constant to third order in t/U, where H 
is given by equation ^ with J2 = 2t'^/U, J^ = 2t'^/U, 
and J3 = 6t'^t'/U^. This is not by chance: the SU{2) 
invariance of -ffpH is automatically inherited to iJeff and 
the chirality is built into the model through the complex 
hopping amplitudes. 

The phase diagram of H on the torus based on Chern 
number computations for a 4 x 5 lattice is depicted in 
figure [5j As the 4x5 lattice is too small to display a 
clear two-fold ground state degeneracy, it may happen 
that the two lowest energy states are not the states that 
resemble the CFT states. We circumvent this problem 
by noting that the Hamiltonian is block diagonal, which 
allows us to select the right states based on quantum 
numbers as stated in the caption. The diagram displays 
the freedom we have to vary the coupling strengths in the 
Hamiltonian while staying within the topological phase. 
In particular, it reveals a ridge with high overlap between 
the CFT state and the ground state of H that connects 
the Hamiltonian studied in the previous sections with the 
region where t,t' <^ U. Choosing, e.g., (j)i = 0.075 x 27r 
and (j)2 — 0.03 x 27r (the white square in the figure), we 
get t/U = 0.10, t'/U = 0.07, and an overlap of 0.985. 

The observation that H and iJpH are closely related 
provides some insight into how our model compares to 
previous proposals. Specifically, Hpa can be seen as a 
sum of two free fermion models, i?kin,t and iJkin,^, and 
a local interaction term. Each of the free models has 
two bands of which the lowest is filled and the highest is 
empty. The Chern number of the filled band (which we 




FIG. 5. Phase diagram. Phase diagram of the Hamilto- 
nian in equation (pi on the torus with J2 = cos((;/>i) cos((;/>2), 
J2 = sin((/!>i) cos(02), and J3 = sin(</>2). C is the total Chern 
number of the states i/'to and i/^xi on a 4 x 5 lattice, where 
^To (V'ti) is the lowest energy state in the subspace spanned 
by all states with the same eigenvalues of S'fot and the trans- 
lation operators in the x- and y-directions as ipTo'^ {4'ti'^)- 
Within the topological phase {C — 1), the two states are well 
separated from higher energy states in the same subspaces 
and flow into each other under flux insertion (like in figure 
[2J3). The background colour gives the overlap between the 
CFT state ([I]) and the ground state of S for a 4 x 5 lattice 
with open boundary conditions. The white circle marks the 
Hamiltonian studied in the first part of the article and the 
white triangle and the white square mark possible parameter 
choices discussed in the text. We omit 4>2 = because the ad- 
ditional symmetries present for this case may cause the lowest 
energy states in the considered subspaces to be degenerate. 



compute using equation (12) in |45] 1 is plus or minus one 
for all nonzero t and t', and the band flattening, i.e., the 
ratio of the gap between the bands to the width of one of 
the bands, is moderate F < 2/(\/2— 1) « 4.8. Our results 
thus show that two copies of a fermionic Chern insulator 
with integer band filling plus local interactions can give 
rise to a bosonic FQH state. This behaviour is different 
from flat band models, where a fermionic (bosonic) Chern 
insulator with flattened and fractionally filled bands plus 
interactions give rise to a fermionic (bosonic) FQH state 
i20k i22i i23i . We therefore conclude that our model pro- 
vides a different mechanism to get the FQH effect in lat- 
tices. 



IV. SIMULATION SCHEME 

The Fermi-Hubbard like model can be simulated in 
fermions in optical lattices as follows. We encode the 
spin up and down states in four internal hyperfinc levels 
as shown in figure [TJi. By choosing the internal states 
appropriately and using polarized light, it is possible to 
construct the optical lattice in figure [!};, where the lev- 
els shown in blue (red) are trapped in a potential with 
minima located at the white (black) squares of a checker- 



board. As detailed in the supplementary information, 
this can be done with eight laser beams travelling in the 
directions (±l,±l,±/3), where /3 is adjustable. 

In this setup, the atoms can tunnel through the po- 
tential barrier between nearest neighbour sites in the 
blue/red lattice, which provides a contribution to the 
next-nearest neighbour hopping terms in the Hamilto- 
nian. We propose to implement the nearest neighbour 
hopping terms by use of laser assisted tunnelling [15] . 
This can be done with the three standing wave laser 
fields in figure [T]d, which we refer to as rl, r2, and b3, 
respectively. As explained further in the supplementary 
information, rl and b3 (r2 and b3) induce nearest neigh- 
bour hops in the horizontal (vertical) direction. A hop 
is accompanied by a Raman transition between internal 
states as illustrated in figure [TJi. The alternating signs of 
the hopping amplitudes originate from the spatial oscil- 
lations of the amplitudes of rl and r2 with a period that 
is twice the lattice spacing, and the j's on the hopping 
amplitudes in the vertical direction appear because the 
amplitude of r2 is imaginary. Note that it is possible to 
choose the internal states in such a way that the energy 
difference between the blue and red spin up states is the 
same as the energy difference between the blue and red 
spin down states, rl, r2, and b3 therefore accomplish the 
hops for both spin up and spin down. 

A particularly convenient feature of the above model 
and implementation scheme is that only on-site interac- 
tions are needed. These occur naturally in the optical 
lattice, and so the Hamiltonian in equation ([3| is im- 
plemented directly. The Hamiltonian ([2| captures the 
physics of ([S]) when t,t' ^ U and the temperature is 
small compared to the gap to the first excited state, 
which is of order t'^/U. The requirements regarding tem- 



perature, interactions, and hopping amplitudes are thus 
the same as those needed to observe the Neel antiferro- 
magnetic order in the standard Fermi-Hubbard model, 
which several groups are pursuing at the moment. Even 
if ([3]) does not exactly reproduce ([2]), it is possible that 
similar physics can be observed. This suggests that, e.g., 
the choice t/U = 1/2 and t' /U — 1/6 (corresponding 
to the white triangle in figure pi), which leads to more 
moderate temperature requirements, may already give 
interesting results. 



ACKNOWLEDGEMENTS 

The authors acknowledge discussions with Fabio Mez- 
zacapo, Hong-Hao Tu, and Shuo Yang. This work has 
been supported by the EU project AQUTE, FIS2012- 
33642, QUITEMAD, and the Severo Ochoa Program. 



AUTHOR CONTRIBUTIONS 

All authors planned the project together and con- 
tributed to the ideas behind the work, AEBN and GS did 
the analytical computations related to the CFT states, 
AEBN did all other analytical and numerical computa- 
tions and drafted the manuscript, AEBN and JIG devel- 
oped the simulation protocol, and all authors contributed 
to the final version of the manuscript. 



COMPETING FINANCIAL INTERESTS 



The authors declare no competing financial interests. 



[1] Tsui, D. C, Stormer, H. L. & Gossard, A. C. Two- [9] 

dimensional magnetotransport in the extreme quantum 
limit. Phys. Rev. Lett. 48, 1559-1562 (1982). 

[2] Yoshioka, D. The Quantum Hall Effect (Springer- Verlag, [10] 
Berlin Heidelberg, 2002). 

[3] Stern, A. Anyons and the quantum Hall effect - A peda- 
gogical review. Annals of Physics 323, 204-249 (2008). [11] 

[4] Stormer, H. L., Tsui, D. C. & Gossard, A. C. The frac- 
tional quantum Hall effect. Rev. Mod. Phys. 71, S298- 
S305 (1999). [12] 

[5] Laughlin, R. B. Anomalous quantum Hall effect: an in- 
compressible quantum fluid with fractionally charged ex- 
citations. Phys. Rev. Lett. 50, 1395-1398 (1983). [13] 

[6] Haldane, F. D. M. Model for a quantum Hall effect with- 
out Landau levels: condensed-matter realization of the 
"parity anomaly". Phys. Rev. Lett. 61, 2015-2018 (1988). [14] 

[7] Wen, X.-G., Wilczek, F. & Zee, A. Chiral spin states and 
superconductivity. Phys. Rev. B 39, 11413-11423 (1989). 

[8] Schroeter, D. F., Kapit, E., Thomale, R. & Greiter, M. [15] 
Spin Hamiltonian for which the chiral spin liquid is the 
exact ground state. Phys. Rev. Lett. 99, 097202 (2007). 



Greiter, M. & Thomale, R. Non-abelian statistics in a 
quantum antiferromagnet. Phys. Rev. Lett. 102, 207203 
(2009). 

Thomale, R., Kapit, E., Schroeter, D. F. & Greiter, M. 
Parent Hamiltonian for the chiral spin liquid. Phys. Rev. 
B 80, 104406 (2009). 

Kapit, E. & Mueller, E. Exact parent Hamiltonian for the 
quantum Hall states in a lattice. Phys. Rev. Lett. 105, 
215303 (2010). 

Nielsen, A. E. B., Cirac, J. L & Sierra, G. Laughlin spin- 
liquid states on lattices obtained from conformal field 
theory. Phys. Rev. Lett. 108, 257206 (2012). 
Greiter M., Schroeter D. F. & Thomale R. Parent Hamil- 
tonian for the non-Abelian chiral spin liquid. Preprint at 
http://arxiv.org/abs/1201.5312vl (2012). 
S0rensen, A. S., Demler, E. & Lukin, M. D. Fractional 
quantum Hall states of atoms in optical lattices. Phys. 
Rev. Lett. 94, 086803 (2005). 

Hafezi, M., S0rensen, A. S., Demler, E. & Lukin, M. D. 
Fractional quantum Hall effect in optical lattices. Phys. 
Rev. A 76, 023613 (2007). 



[16] Bloch, I., Dalibard, J. & Zwerger, W. Many-body physics [35 

with ultracold gases. Rev. Mod. Phys. 80, 885-964 

(2008). 
[17] Cooper, N. R. Rapidly rotating atomic gases. Advances [36 

m Physics 57, 539-616 (2008). 
[18] Tang, E., Mei, J.-W. & Wen, X.-G. High-temperature 

fractional quantum Hall states. Phys. Rev. Lett. 106, 

236802 (2011). [37 
[19] Sun, K., Gu, Z., Katsura, H. & Das Sarma, S. Nearly 

flatbands with nontrivial topology. Phys. Rev. Lett. 106, 

236803 (2011). [38 
[20] Neupert, T., Santos, L., Chamon, C. & Mudry, C. Frac- 
tional quantum Hall states at zero magnetic field. Phys. 

Rev. Lett. 106, 236804 (2011). [39 

[21] Roy, R. & Sondhi, S. L. Fractional quantum Hall effect 
without Landau levels. Physics 4, 46 (2011). 

[22] Sheng, D. N., Gu, Z.-C., Sun, K. & Sheng, L. Fractional [40 
quantum Hall effect in the absence of Landau levels. Nat. 
Commun. 2, 389 (2011). [41 

[23] Wang, Y.-F., Gu, Z.-C., Gong, C.-D. & Sheng, D. N. 
Fractional quantum Hall effect of hard-core bosons in 
topological flat bands. Phys. Rev. Lett. 107, 146803 [42 
(2011). 

[24] Regnault, N. and Bernevig, B. A. Fractional Chern Insu- 
lator. Phys. Rev. X 1, 021014 (2011). [43 

[25] Nayak, C, Simon, S. H., Stern, A., Freedman, M. & Das 
Sarma, S. Non-Abelian anyons and topological quantum 
computation. Rev. Mod. Phys. 80, 1083-1159 (2008). 

[26] Parameswaran, S. A., Roy, R. & Sondhi, S. L. Fractional [44 
Chern insulators and the Woo algebra. Phys. Rev. B 85, 
241308(R) (2012). 

[27] Yao, N. Y., et al. Realizing fractional Chern [45 
insulators with dipolar spins. Preprint at 
|http://arxiv.org/abs/12l"2:483 9 (2012). [46 

[28] Kalmeyer, V. & Laughlin, R. B. Equivalence of the 
resonating-valence-bond and fractional quantum Hall 
states. Phys. Rev. Lett. 59, 2095-2098 (1987). [47 

[29] Laughlin, R. B. Spin hamiltonian for which quantum hall 
wavefunction is exact. Annals of Physics 191, 163-202 
(1989). [ 

[30] Bloch, I., Dalibard, J. & Nascimbene, S. Quantum sim- 
ulations with ultracold quantum gases. Nature Phys. 8, 
267-276 (2012). 

[31] Jordens, R., Strohmaier, N., Giinter, K., Moritz, H. & [49 
Esslinger, T. A Mott insulator of fermionic atoms in an 
optical lattice. Nature 455, 204-207 (2008). [50 

[32] Schneider, U. et al. Metallic and insulating phases of re- 
pulsively interacting fermions in a 3D optical lattice. Sci- 
ence 322, 1520-1525 (2008). [51 

[33] Trotzky, S. et al. Time-resolved observation and control 
of superexchange interactions with ultracold atoms in 
optical lattices. Science 319, 295-299 (2008). [52 

[34] Nielsen, A. E. B., Cirac, J. I. & Sierra, G. Quantum 
spin Hamiltonians for the SU{2)k WZW model. J. Stat. 
Mech. P11014 (2011). 



Moore, G. & Read, N. Nonabelions in the fractional 
quantum Hall effect. Nuclear Physics _B360, 362-396 
(1991). 

Thouless, D. J., Kohmoto, M., Nightingale, M. P. & 
den Nijs, M. Quantized Hall conductance in a two- 
dimensional periodic potential. Phys. Rev. Lett. 49, 405- 
408 (1982). 

Niu, Q., Thouless, D. J. & Wu, Y.-S. Quantized Hall 
conductance as a topological invariant. Phys. Rev. B 31, 
3372-3377 (1985). 

Kohmoto, M. Topological invariant and the quantization 
of the Hall conductance. Annals of physics 160, 343-354 
(1985). 

Hatsugai, Y. Characterization of topological insulators: 
Chern numbers for ground state multiplet, J. Phys. Soc. 
Jpn. 74, 1374-1377 (2005). 

Thouless, D. J. Level crossing and the fractional quantum 
Hall effect. Phys. Rev. B 40, 12034-12036 (1989). 
Lauchli, A. Numerical Simulations of Frustrated Sys- 
tems. Springer Series tn Solid-State Sciences 164, 481- 
511 (2011). 

Jiang, H.-C, Wang, Z. & Balents, L. Identifying topo- 
logical order by entanglement entropy. Nature Phys. 8, 
902-905 (2012). 

Fendley, P., Fisher, M. P. A. & Nayak, C. Topologi- 
cal entanglement entropy from the holographic partition 
function. .Journal of Statistical Physics 126, 1111-1144 
(2007). 

Bravyi, S., DiVincenzo, D. P. & Loss, D. Schrieffer- Wolff 
transformation for quantum many-body systems. Annals 
of Physics 326, 2793-2826 (2011). 

Shi, T. & Cirac J. I. Topological phenomena in trapped- 
ion systems. Phys. Rev. A 87, 013606 (2013). 
Jaksch D. & Zoller P. Creation of effective magnetic fields 
in optical lattices: the Hofstadter butterfly for cold neu- 
tral atoms. New J. Phys. 5, 56 (2003). 
Haldane, F. D. M. Many-particle translational symme- 
tries of two-dimensional electrons at rational Landau- 
level filling. Phys. Rev. Lett. 55, 2095-2098 (1985). 
Hastings, M. B., Gonzalez, I., Kallin, A. B. & Melko, R. 
G. Measuring Renyi entanglement entropy in quantum 
Monte Carlo simulations. Phys. Rev. Lett. 104, 157201 
(2010). 

Kitaev, A. & Preskill, J. Topological entanglement en- 
tropy. Phys. Rev. Lett. 96, 110404 (2006). 
Levin, M. & Wen, X.-G. Detecting topological order in a 
ground state wave function. Phys. Rev. Lett. 96, 110405 
(2006). 

Haque, M., Zozulya, O. & Schoutens, K. Entanglement 
entropy in fermionic Laughlin states. Phys. Rev. Lett. 
98, 060401 (2007). 

Zhang, Y., Grover, T. & Vishwanath, A. Topological en- 
tanglement entropy of Z2 spin liquids and lattice Laugh- 
lin states. Phys. Rev. B 84, 075128 (2011). 



Supplementary information 

CFT STATES FROM THE SU{2)i WESS-ZUMINO-WITTEN MODEL 

Background 

The Wess-Zumino-Witten (WZW) model based on the Kac-Moody group SU{2)i is a CFT with central charge 
c = 1 and two primary fields 4>o{z) and 4>i/2iz) with conformal weights /iq = and /ii/2 = 1/4, respectively. The 
field 4>j{z) is associated with the total spin j and has components ipj^m^z) with m ~ — j, — j + 1, . . . , j. From these 
fields one can construct chiral correlators, also called conformal blocks, 

and interpret the result as a quantum state IV'fe^"'") = X^m m-, m V'fe^"'"('TT'ii '7^2, . . . , rrijv) Imi, ^2, . . . , niM) of 
TV particles with spins j„ at fixed positions Zn- (•) stands for vacuum expectation value, the word chiral refers 



to the fact that only z„ and not z„ appears on the right hand side of (SI), and the index k takes into account 
that there may be more than one conformal block for a given choice of fields. The number of conformal blocks 
on a given Riemann surface with genus g is dictated by the fusion rules of the primary fields. Specifically, for N 
primary fields (j)i/2 (which is the interesting case since ^o represents the identity) the number of conformal blocks 
is 2^ if N is even, whereas the chiral correlator vanishes if N is odd. Hence, on the sphere, i.e. g = 0, there is only 
one chiral correlator, while on the torus, i.e. 5 = 1, there are two. 

CFT state in the plane 

The plane geometry considered in the main text is the Riemann sphere. The chiral correlator is well-known in 
this geometry, so here we simply note that the CFT expression leading to equation (1) in the main text is 

V'po'^(si,S2,--- ,Sn) OC (01 iyL(2;i)0i £a(z2) .. .01 £|l(zjv)) (S2) 

with N even. 

CFT state on the cylinder with closed ends 

The plane (or more precisely Riemann sphere) can be mapped to the cylinder with closed ends and circumference 
Ly via the coordinate transformation z — g^'^(h+^h)/Ly ^ where z is a point in the plane written as a complex 
number and (IxJy) is the corresponding point on the cylinder (note that ly is periodic with period Ly, whereas 
Ir^ is not). The cylinder has closed ends, and hence genus zero, because e^'^^'^+^'s')/^'' is a single point rather 
than a circle for l^ = ±oo. The wavefunction V'co"^ '^^ ^^'^ cylinder is therefore easily obtained from V'po""'" 
by choosing z„ = e'^'^i'-^.'^+'-'-v.'iJ/'^v for n — 1,2, ... ,N. For the particular case of a square lattice, we choose 
l.,n G {-{L. - l)/2, -{L, - l)/2 + 1, . . . , (L, - l)/2} and ly^r. G {1, 2, . . . , L J. 

CFT states on the torus 

To derive the conformal blocks on the torus, we use the fact that the SU{2)i WZW model can be bosonized in 
terms of a massless free scalar field (f{z,z) compactified on a circle of radius R = ^2. This leads to the result 



n 



I k 



k 



(S3) 
where '^{zi) {Lp{zi)) is the chiral (antichiral) component of the scalar field ^{z, z) — (p{z) + (p{z) and : • : denotes 
normal ordering. An expression for the correlator on the left hand side has been obtained in [Ij , and by rewriting 



this expression we can extract the conformal blocks ( SI I up to phase factors that are independent oi Zi, Z2, ■ ■ ■ , Zi\[. 
The phase factors are determined by noting that the conformal blocks on the torus reduce to the conformal block 
on the plane if all the spins sit within a small region. 

Before writing down the result for the conformal blocks, we note that a torus can be defined by specifying two 
complex numbers Wi and W2 and identifying points in the complex plane that differ by nuji + muj2, where n and 
m are integers. (The appropriate choice for the square lattice, e.g., is uji = L^ and uj2 — iLy.) The modular 



parameter r = ^2/^1 is assumed to have positive imaginary part without loss of generahty. Defining the scaled 
positions w„ = z„/a;i, the conformal blocks read 



^r^(.i, S2, ...,s^)^ ^(_i)5::r..(«-i)(..+i), - ,,^„ 







2t\ H E{wn,w^y-'-^'\ (S4) 



where k takes the values and 1/2, 



1 for 2„s„ 
otherwise 



(S5) 



r]{q) is the Dedekin eta function, 



00 



^(9) -9'/'' 11(1 -9")' g = e2--, (S6) 

n=l 

(^a,biw,T) is the theta function, 

00 



E(wn,Wm) is the prime form on the torus. 



9l l(Wn — Wm,T) d0 1 i(w,t) 



2 ' 2 



9\ 1(0, t) ' 2^2^ ' ' rfw 



(S8) 



10=0 



and the operation (•)''n'''"/2 jg jjj^fjg single valued by always choosing the phase of E{wn,Wm) to be within the 
interval ] — tt, tt]. 

Let us consider the special case of an L^ x Ly square lattice and choose the origin such that J2n ^n ~ ^■ 
If we regard a spin up as an occupied site and a spin down as an empty site, we find after several analytical 



manipulations that ( S4 1 is proportional to the corresponding continuum FQH states provided in [21 [3] except that 
the possible positions of the particles are restricted to the lattice and except for some sign factors that make the 
state into a spin singlet. In other words, the CFT states are precisely the Kalmeyer-Laughlin states on the torus. 
We note also that the states in [3] allow for twisted boundary conditions and therefore provide a natural way to 
generalize the CFT states given below to this case. 



The CFT states on the torus are obtained as linear combinations of the conformal blocks ( S4 ) chosen in such a 
way that the CFT states refiect the symmetries of the Hamiltonian in equation (2) in the main text. Specifically, 
we would like the CFT states to be eigenstates of the translation operators, T^ and Ty, in the x- and y-directions. 
Acting with T^ and Ty on the conformal blocks, we find after a lengthy computation that the appropriate choice 
for an even-by-odd lattice is 

'4>^q'^{si,S2, . . . , Sn) OC ■0^^'^(s1, S2, . . . , S n) 

(even- by-odd lattice). (S9) 

'4'Tl'^isi,S2, . . . , Sn) OC 'lJj'^^^{si,S2, . . . , Sn) 

For the even-by-even lattice with L^ ~ Ly both of the conformal blocks are invariant under translations, so we 
need to consider additional symmetries. Requiring the CFT states to be eigenstates of the operator that rotates 
the lattice by ninety degrees, we get 

^To'^isi,S2,-. .,SAr) OC V^f ^"^(si , S2, . . . , Sjv) + (-1 + V^) Ip'^J'^ (si, S2, ■ ■ ■ , Sn) 
V'Tf'^(si,S2,-- -^Sn) OC V'o'^'^(si,S2,...,SAr) + (-1 - V^} tp'^J'^ {si, S2, ■ ■ ■ , Sn) 

This construction also guarantees the CFT states to be orthogonal, which is not necessarily the case for the 
conformal blocks. 



10 
OPTICAL LATTICE POTENTIAL 

The main idea behind the scheme we propose to reahze the optical lattice potential in figure Ic in the main text 
is to encode the blue (red) states in internal levels that interact more strongly with right (left) circularly polarized 
light than with left (right) circularly polarized light and then create a checkerboard pattern of regions with highest 
intensity of right circularly polarized light and regions with highest intensity of left circularly polarized light. Since 
the blue and red sublattices are square lattices when viewed from the {x + y)- or (x — j/)-directions, it is natural to 
let laser beams travel along these directions to create a standing wave pattern of intensity minima and maxima. 
In order to be able to control the strengths and phases of right and left circularly polarized light in the beams 
independently, however, the momentum vectors of the beams should have a nonzero z-component, which leads to 
the eight travelling directions (±l,±l,±/3) mentioned in the main text. We note that /3^ « 2 is a particularly 
convenient choice because the trapping lasers then have approximately the same frequency as the fields used for 
inducing hops and one can use the same set of excited states for implementing the trap and for implementing the 
hopping terms (the frequencies should, of course, not be exactly the same to avoid interference). 

Let us demonstrate explicitly how this can be done for fermions with one valence electron, when the ground state 
manifold is a ^5*1/2 orbital, and the light fields couple the ground states off-resonantly to a ^Pi/2 orbital (see figure 



Sli), which can be achieved with alkali atoms. We first create a standing wave pattern in the (x + j/)-direction 



by adding up the fields 



Eia{r,T) = Re (e l^/2e+ + a\/2ze_ + ^(1 + i){l - a)e\ e^^A'-'^+y^f^k^-^-^'^iA ^ (g^^) 



Eib{r,T) = Re (e (V2e+ - aV2ie- - ^(1 + i){l + a)eA 



— ik^{x+y) ik^z —iuJiT 



e • "' e e 



(S12) 



where r = {x,y, z), T is time, E — E* , e+ — (— 1, — z, 0)/\/2 is the polarization vector of right circularly polarized 
light, e_ = (1, —i, 0)/-\/2 is the polarization vector of left circularly polarized light, Sz = (0, 0, 1) is the polarization 
vector of z-polarized light, a is a real parameter, 

' ' '^ fc=^, (S13) 



^ y2T^' ' y2T^' A 

A is the wavelength of the fields, and tui is the frequency. Note that the intensity profile of the right (left) circularly 
polarized component is proportional to cos'^{kx{x + y)) (sin {kx{x + y))) and that the coeflficients of £*z have been 
chosen such that the polarizations of the fields are perpendicular to the directions of travelling as it should be for 
light fields. As seen in the figure, the presence of the z-polarized component causes the trapping lasers to induce 
Raman transitions between the two ground states. This undesired effect can be cancelled by adding the fields 

E2a{r,T) = Re (e (V2e+ + aV2ie- - ^(1 + «)(1 - a)e,') e^fc^("+^)e-*'=-"e-*"=^') , (S14) 

E2t{r,T) = Re (e (V2e+ - aV2ie^ + ^(1 + «)(1 + a)e,') e-«fex(:>.+y)g-^fc.^g-^<^.^^ ^ ^g^^^ 

where the frequency uj2 is assumed to be slightly different from uji such that the fields do not interfere coherently. 
The fields 

E3a{r,T) = Re (e (V2e+ + aV2te^ - ^{1 - i){l + a)ez^ e-'''-i^^-y)e^k^^^-^-',T\ ^ (g^g^ 

E3b{r,T) = Re (e (V2e+ - aV2ie- + ^(1 - i)(l - a)e,^ e^'^-^^^^^e^'^-^e-^'^^^') , (S17) 

E4a{r,T) = Re (e (V2e+ + aV2ie- + ^(1 - i){l + a)e^ j e-»'=.(^-!/)e-*'=-^e-*'^^^ j , (S18) 



Eibir,T) = Re (e (V2e+ - aV2ie. - ^(1 - i){l - a)s^ e'''-'^^-y'>e-'''^'e~'^^A , 



(S19) 



do the same for the {x — y)-direction, where W3 and lu^ are slightly different from uji and uj2 to avoid coherent 
interference. 



■5 = 1/2 |e_) |e+) 

L = l -r 

J = 1/2 I A 




S = 1/2 
L = 

^ = 1/2 \ |g_) \g+) 



It) lt> 



i;> 



' 2 ' I 2' 2' 



11 



|3 -3\ |3 -1\ |3 1\ |3 3 



/ I 2' 2/ I 2 ' 2/ 



FIG. SI. Atomic levels and light-atom interaction, a, Ofl-resonant coupling of a S1/2 ground state manifold to a 
^Pi/2 excited state manifold with polarized light. S is the spin angular momentum, L is the orbital angular momentum, 
and J is the total angular momentum obtained by coupling S and L. b, Hyperfine levels of the ground state manifold when 
the nuclear spin is / = 1. The states are labelled \F,mF), where F is the total angular momentum obtained by coupling 
J and 7, and ttif is the projection of the total angular momentum on the z-axis. The choice of blue and red spin up and 
down states is indicated. 



Eliminating the excited states with a Schrieffer- Wolff transformation, we get the trapping potentials 

V- - -Vo{2l3^ + a^- l)(sin2(A:,(a; + y)) + sm^ik,{x - y))) + 2^0(2/3^ + a'), (S20) 



V+ = -V^o(l + 2a2^2 _ a^)(^cos\k,{x + y)) + cos\k^{x - y))) + 21/o(l + ^a^P^), 



(S21) 



for 



and 



(?+), respectively, where Vq oc E /{f3 A) and the detuning A, defined in the figure, is negative 
(positive) for red (blue) detuning. For — Vb(2/3^ + q^ - 1) > and -Vb(l + 2a^/3^ — a^) > 0, we thus conclude 
that \g-) is trapped around the minima of (sin {kx{x + y)) +sin {kx{x --y))), whereas \g+) is trapped around the 
minima of (cos^(fc2;(a; + y)) + cos^{kx{x — y))), which gives the desired checkerboard pattern. 

So far, we have considered the states at the level of fine structure. Assuming, as an example, that the nuclear 
spin is / = 1, which is the case for ^Li, the ground state manifold in fact consists of six states as depicted in figure 
|S1|3. The trapping potentials of these states are 





Va 
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amounts to 





-V^o(l + 2a2/32 - a'){cos\k,{x + y)) + cos2(fc,(x - y))) + 21/o(l + '^a'^'^), 

-K) Qa2/32 -a' + l- "^pA (cos2(fc,(x + y)) + co^\kx{x - y))) + 2V^ (l + \(^^P''^ 

-Fo (o? - ^a2/32 + ^/32 _ l") (sin2(fc,(a; + y)) + sm\k,{x - y))) + 2V^ L? + ^p'\ , 
-K)(a' + 2/32 - l){s^n\kx{x + y)) + sin2(fc,(a: - y))) + 2Vo{o? + 2/3^). 



(S22) 
, (S23) 

(S24) 
(S25) 



+ 1 - I/32) > and -Vq (a^ - fa^^^ ^ |^2 _ ^-^ > q^ .^j^j^^j^ f^j. ^.gj detuning {Vq < 0) 



/32 



1 



< a < 



P2 



(or 1/5 < a^ < 5 for p' = 2), 



(S26) 



we can get the desired potential by choosing the blue and red spin up and down states as shown in figure pTp. 
Since the potentials shift the two blue states by the same amount and the two red states by the same amount, 
the energy difference between the blue and red spin up states is the same as the energy difference between the 
blue and red spin down states, which is what allows us to drive the hops of spin up and spin down with the same 
lasers. If this symmetry cannot be achieved in a given setup, however, one can double the number of standing 
wave laser fields used for implementing the hops and drive the transitions independently. The energies of the 
states ||, |) and ||, — f ) are, in general, shifted by different amounts so that transitions to these states induced by 
the hopping lasers can be avoided automatically. Note also that the difference between the maxima and minima 
of the potentials can be adjusted. For (3' = 2, e.g., the difference is | — 2Vb(— |a^ + 1)1 ^i' the blue lattice 



potential and | — 2Vb(|a^ "1)1 ^'^^' the red lattice potential. Therefore a^ and /3^ allow us to control the ratio 
between the tunnelling rates between nearest-neighbour sites in the blue and red lattices. Finally, the energy of 
a fermion sitting on a blue lattice site must be different from the energy of a fermion sitting on a red lattice site 
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to make the proposed implementation of the hopping terms work. If the contributions to this difference coming 
from the optical lattice potential and the field labelled b3 in the main text are not appropriate, one can control 
the difference by adding an additional laser field along the 2;-axis with right or left circular polarization. 

LASER INDUCED HOPPING TERMS 

The transitions induced by the laser fields rl, r2, and b3 in the main text can be understood by noting that the 
transition amplitudes are proportional to the spatial integral of the product of the Wannier functions at the two 
sites and the two field components that are involved in the transition. Since the Wannier functions decay rapidly, 
there are only transitions between nearby sites. 

Let us first consider the jumps between neighbouring blue and red lattice sites induced by rl and b3. The 
transition amplitude for jumps along the y-axis vanishes because the Wannier functions and bS are even functions 
and rl is an odd function with respect to reverting the x-axis around the x-coordinate of the sites. Jumps along the 
X-axis are allowed, and since the sign of the amplitude of rl between two neighbouring sites is alternatingly plus 
and minus, the signs of the hopping amplitudes also alternate. r2 and b3 similarly induce hops in the y-direction, 
but not in the x-direction. For the choice of hyperfine levels considered above, it is the right (left) circularly 
polarized component of bS that is involved when a fermion in the spin up (down) state jumps. Therefore one can 
make the hopping amplitudes for up and down spins equal by adjusting E± appropriately, i.e., by adjusting the 
polarization state of b3. 

Let us also consider transitions involving one photon from rl and one photon from r2. In this case, the fermion 
cannot change its internal state due to energy conservation, and it cannot jump to a site of different colour. If 
the fermion remains at the same site, the transition amplitude is zero due to symmetry. If the transition involves 
a jump to a next-nearest neighbour site, there is destructive interference between absorbing a photon from rl 
and emitting a photon into r2 or absorbing a photon from r2 and emitting a photon into rl, and therefore no 
transitions occur. 

rl, r2, and bS also induce transitions, where a fermion absorbs and emits a photon from the same beam without 
changing its internal state and either remains at the same site or jumps to a next-nearest neighbour site. The 
former case gives rise to the trapping in the z-direction and also slightly modifies the potential in the xy-plane. 
The modification can, however, be made insignificant by decreasing E and increasing E± without affecting the 
nearest-neighbour hopping amplitudes that are proportional to EE± . The latter case gives rise to a contribution 
to the next-nearest neighbour hopping terms in the Hamiltonian. For a given choice of E± it may happen that 
hops from blue to blue lattice sites occur at a different rate than hops from red to red lattice sites. As mentioned 
above, however, one can lower the potential barriers between the sites for the sublattice with slow jumps relative 
to the other to compensate. 
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